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Abstract. We prove that the first eigenvalue of the Dirichlet Laplacian 

2 j 2 

for a triangle in the plane is bounded above by w gAi , where L is the 
perimeter and A is the area of this triangle. We show that the constant 9 
is optimal and that the optimal constant for the lower bound of the same 
form is 16. This gives a positive answer to a conjecture made in 0. 



The purpose of this paper is to prove the following theorem. 

Theorem 1.1. Let T be a triangle in a plane of area A and perimeter L. 
Then the first eigenvalue Ay of the Dirichlet Laplacian on T satisfies 



The constants 9 and 16 are optimal. 

The lower bound was proved in a more general context in [31- In Section 
6 we show that for "tall" isosceles triangles there is an asymptotic equality 
in the lower bound. Hence it is impossible to decrease the constant 16. 

The upper bound was recently stated as a conjecture in |2j and numerical 
evidence for its validity are given in pQ. Bounds of this form but with 
different constants have been the subject of many papers in the literature. 
The eigenvalue of any doubly-connected domain is bounded above by the 
same fraction but with the constant 4, see [§] and remarks in 0. There is 
also a sharper upper bound due to Freitas ( t 2 J which is not of this form but 
it seems that in the worst case ( "tall" isosceles triangle) it gives the constant 
6, in the best (equilateral) 9. It is worth noting that the constant 9 may not 
be improved since equilateral triangles give equality in the upper bound of 
the Theorem ll.il 

The spectral properties of a Dirichlet Laplacian on an arbitrary planar 
domain are important both in physics and in mathematics. Unfortunately, 
it is almost impossible to find the exact spectrum even for some simple 
classes of domains. Except for rectangles, balls and annuli, not much can 
be said in general. In the case of triangles the full spectrum is known only 
for equilateral and right triangles with smallest angles tt/4 or tt/6. For more 
information about these we refer the reader to For all other triangles, 
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the best we can hope to do is to give bounds for the eigenvalues, such as 
those given above. 

Even though Theorem 11.11 gives sharp bounds in the sense that the con- 
stants are the best possible given the form of the bound, there is certainly 
room for improvements. In fact, sharper lower bounds are already known, 
see j2j. One of these bounds is good for both equilateral and "tall" trian- 
gles. It gives the constant 9 for the first and 16 for the second. The upper 
bound good in both cases is still unknown. To the best of our knowledge 
it is also not clear what is the correct bound for the isosceles triangle with 
base almost equal to the half of the perimeter, but we think it should be 16. 

By comparing our numerical results with the numerical studies contained 
in pP, Section 5.1, we conjecture that 

Conjecture 1.2. Let T be a triangle in a plane of area A and perimeter L. 
Then the first eigenvalue Xt of the Dirichlet Laplacian on T satisfies 

. . 7T 2 L 2 7V37T 2 7T 2 L 2 V3TT 2 

(L2) 16^ + ^2l"- AT - V2A 2+ ^A- 

Here both bounds are of the form 

considered in 

The lower bound, with 9 = 7r 2 /16, is the best bound we can expect 
given this particular form. Indeed, this is the only bound which is sharper 
then the lower bound of Theorem 11.11 and which might be true for "tall" 
triangles. The upper bound from our main result is also of this form, but 
with 9 = 7r 2 /9. Hence the conjectured upper bound is sharper (9 = 7r 2 /12), 
and it is the best in the sense that the bound with 9 = 7r 2 /13 is not valid. 
But, since only the constant 16 can give a good upper bound for "tall" 
triangles, it is not possible to find a bound of the form £"3 which is good for 
both equilateral and "tall" triangles. 

Our proof of the upper bound from the Theorem 11.11 contains two main 
parts. The first deals with an "almost equilateral" triangles. That is, with 
the triangles for which the longest side is comparable to the shortest side. 
For these our strategy is to find a suitable test function ip. That is, we try 
to find a function which is on the boundary of the triangle T and apply 
the Rayleigh quotient to get the upper bound for Xt- We get 

(1-3) A T <_ 

Jt ^ 

This part of the proof is contained in Sections 2 to 5. Included in Section 2 
are also some preliminary results. 

The second part of the proof, contained in Section 6, deals with "tall" 
triangles. These can be approximated by a circular sections for which the 
eigenvalues can be found explicitly. 
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2. ElGENFUNCTIONS AND NOTATION 



An arbitrary triangle T' can be rotated and rescaled to obtain a triangle 
T with vertices (0, 0), (1,0) and (a, b). This, together with the fact that the 
bound in the main theorem is invariant under translations, rotations and 
scaling allow us to restrict our attention to the triangles with such vertices. 
We can also assume that the side contained in the rr-axis is the shortest. 
Hence we have that 

a 2 + b 2 > 1 and a < 1/2 

for our triangles. We will denote the length of the other two sides by M and 
N, with N denoting the longest. 

We start with the first eigenfunction of an equilateral triangle, and we 
will proceed as in |2j. Such function is given by 
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We can compose / with a linear transformation to obtain a function 
which is equal to on the boundary of T. Namely consider 
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This function was used in [2] to obtain the upper bound from the Rayleigh 
quotient. Since the function / is the first eigenfunction of the Dirichlet 
Laplacian on an equilateral triangle and its eigenvalue gives equal sign in 
the main bound, it is reasonable to expect that by taking any linear trans- 
formation we can only decrease the constant 9 in the Theorem ll.il 

Hence we want to find another eigenfunction of some other triangle. We 
will use the eigenfunctions of the equilateral triangle to find a test function 
for the right triangle with angles 7r/3 and tt/6. In the recent paper [3] the 
author constructs two families of eigenfunctions of the equilateral triangle. 
The antisymmetric mode has the property that it is on the altitude. Thus, 
such a function is also the eigenfunction for the right triangle. We can 
then take the antisymmetric eigenfunction corresponding to the smallest 
eigenvalue as our test function. A calculation leads to the following function 



(2.3) 



This function, as can be easily checked, is in fact the eigenfunction of the 
Dirichlet Laplacian on the triangle with the vertices (0, 0), (1, 0) and (0, 



g(x, y) = sin ^V37ryJ sin (ttx/3) 

+ sin (^VSiry/S^ sin (57rx/3) 
+ sin (2v / 37r?//3) sin (Attx/3) 
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The corresponding eigenvalue gives the better bound than the one in the 
Theorem ll.il the constant is about 9.6. Therefore, a linear transformation 
of this function should give a correct bound at least for the neighborhood of 
the point (0, \/3). By applying a suitable linear transformation we get the 
second test function 



(2.4) 




+ sm 



Similarly, we can obtain the last two test functions. One will be a linear 
transformation of the eigenfunction of the triangle with vertices (0, 0), (1, 0) 
and (1, v3)- The other a linear transformation of the eigenfunction of the 
triangle with the vertices (0,0), (1,0) and (0, I/a/3)- We get: 
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Now we can take a linear combination of these test functions. That is, 
we consider 



(2.6) 



if>(x,y) = aupi(x,y) + /3<p 2 (x,y) +jip 3 (x,y) +e<f>(x,y). 



and we can calculate the Rayleigh quotient for this function. After optimiz- 
ing over all possible values of a, f3, 7 and e, this will give an appropriate 
bound for the first eigenvalue. To prove Theorem II. II we have to check that 
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for some a, (3, 7 and e (possibly depending on T). The last inequality is 
equivalent to 

(2.8) 9A 2 [ |VVf < vr 2 L 2 f ip 2 . 

Jt Jt 

Since the function ip is given explicitly and is a trigonometric function, 
it is possible to find the exact values of this integrals but calculations are 
very cumbersome. For this reason we will do the long calculations in Mathe- 
matica. However, we wish to emphasize the fact that all the calculations 
are done symbolically. 

By our assumptions we have L = 1 + Va 2 + b 2 + \/ (a — l) 2 + b 2 and 
A = 6/2. As a result of running Mathematica we get that to prove the 
inequality ()2.8j) we have to find a, (5, 7 and e such that the inequality 

(2.9) 

> 8041366333x 

I (-1594323 - 1792090a + 531441(a 2 + b 2 ) + 201600 (3 + a 2 + b 2 ) it 2 ) a 2 

+ (-2854972 + 729208a + 531441(a 2 + b 2 ) + 201600 (3 + (a - l) 2 + 6 2 ) tt 2 ) 1 

+ (531441 - 1792090a - 1594323(a 2 + b 2 ) + 201600 (l + 3a 2 + 36 2 ) tt 2 ) 7 2 | 

+ 5558192409369600 (l - a + (a 2 + 6 2 )) vr 2 e 2 
+ 67672797192x 

I (729V3 (454 - 128a + 339(a 2 + b 2 )) + 24640 (4 - 8a + 9(a 2 + b 2 )) tt) 76 
+ (729v / 3 (665 - 780a + 454(a 2 + b 2 )) + 24640 (5 + 4(a 2 + b 2 )) tt) pe 
+ (729V3 (339 - 128a + 454(a 2 + b 2 )) + 24640 (9 - 8a + 4(a 2 + b 2 )) wj ae\ 

+ (1990033124626008a + 2553294638054160\/3 (3 - 2a + 3(a 2 + b 2 )) tt) 
+ 1151172000 (35341051 - 26756686a + 32479596(a 2 + b 2 )) /3 7 
+ 189 1819452341268271 + -73323642839420a + 73323642839420(a 2 + b 2 ) 
- 7993561 0875 120V37T + 24336134222400^3 (-a + a 2 + b 2 ) 7r|a/3 

- 9 ^1 + \/{-l + a) 2 + b 2 + Va 2 + b 2 ^j |444001222376712v^ (a + 7) e 

- 1629547920tt (VSa (4251/3 - 99484 7 ) - 113696 (a + [3 + 7) e) 
+ 51464744531200tt 2 (a 2 + /3 2 + 7 2 + 2e 2 ) 

+ 3/3 (346474423262177a + 85272 (3297684500 7 + 1735627257V3e)) } 
is valid. 

This expression clearly shows that it would be very difficult to do the 
calculations by hand. Notice that this expression depends only on b 2 and a. 
Also note that the "building" blocks for the expression are exactly equal to 
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the length of the sides of the triangle T. Hence we make the substitution 
M 2 = a 2 + b 2 and N 2 = (a — l) 2 + b 2 . As a result we get a polynomial of 
degree 2 in M and JV, where JV > M > 1. For further simplification (and 
to improve our chances of finding the appropriate coefficients) we divide all 
triangles into 4 classes. Each class will be handled in a separate section 
whose number corresponds to the case: 

3) The triangles with JV > 2 and M < 15, 

4) 1 < JV < 2 and (JV + l)/2 < M < 2, 

5) 1 < JV < 2 and 1 < M < (JV + 1) /2, 

6) M > 15. 

The method used to handle the last case will be totally different than the 
previous ones. 

3. Case: JV > 2 and M < 15 

Let us take e = (3 = 0, a = 1 and 7 = — 1/6. This simplifies (|2.9|) to 
(3.1) 

>P(M, JV) = -90851035780 - 16374894040M 2 + 339299845937V 2 
- 272432160\/3 (10 - 8M + 10M 2 - 8JV - 8MN + 3A^ 2 ) ir 
+ 28828800 (689 - 148M + 199M 2 - 148A^ - 148MiV - 74iV 2 ) vr 2 . 

To show this inequality we first find all the critical points of the right side 
and later check the values on the boundary. Both 8mP and djyP are linear 
with respect to M and N, therefore we have exactly 1 critical point with 
N ~ —42.2. Hence it is enough to check this inequality on the boundary 

The boundary conditions are given by M = JV, M = 15, N = 2 and 
M = N — 1. For each of these P is a quadratic equation and we just have 
to check that the roots are outside of the bounds for M or N and that the 
inequality is true at the endpoints. We have 

• P(M,M) = for M ps 1.6 and M m 15.15; P(2,2) < and 
P(15,15) < 0, 

• P(15,N) = for N ss 14.97 and N w 42.5; P(15, 15) < and 
P(15,16) < 0, 

• P(M, 2) = for M w 0.96 and M ps 2.61; P(l, 2) < and P(2, 2) < 
0, 

• P(N - 1, N) = for JV w 1.97 and JV w 20.56; P(l,2) < and 
P(15, 16) < 0. 

This shows that the desired inequality is true on the boundary and therefore 
everywhere. 

4. Case: 1 < N < 2 and (JV + l)/2 < M < 2 

In the next 2 sections we will have to deal with cases for which an equi- 
lateral triangle (JV = M = 1) is one of the possible triangles. We take 
e = 1, since only the eigenfunction of the equilateral triangle can give the 
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constant 9 in the Theorem 11.11 We also need for all the other coefficients 
to vanish near the equilateral triangle. Let us take 7 = and a = /3. We 
just have to choose the common value for a and (3. Due to the nature of 
the already very complicated calculations we cannot afford to pick a very 
complicated coefficient, thus we take a = (N + M — 2)/2. This choice has 
one additional advantage. In this case we are working with the eigenfunc- 
tions of the following triangles: One equilateral and two right triangles with 
shortest side (0, 0) — (1, 0). Hence we have a symmetry about a = 1/2, or in 
terms of M and N, about M = N. Therefore it is natural to introduce the 
rotated coordinates U = (M + N)/2 - 1 and V = (N - M)/2. Note also 
that a = j3 = U . This also moves the equilateral triangle to the origin. 
After applying these transformations the inequality (|2.9|) becomes 

(4.1) 

>P(U, V) = U 2 (3293385188722144 - 451048860827136^3 - 952832984463360tt 

- 874782993324240V3tt + 463182700780800tt 2 - 4817666363010084/7 

+ 916192998555120V3C/ + 710514087323040v / 3tt[/ - 330844786272000tt 2 [/ 

- 1072431834636645C/ 2 + 346350633108480V37TL/ 2 - 33084478627200tt 2 [/ 2 ) 
+9V^ 2 (441126381696007r 2 + 355514206276944V3J7 + 105870331607040tt?7 

+ 177818984344461C/ 2 + 36504201333600v / 3vr[/ 2 + 25732372265600tt 2 [/ 2 ) 

This is a polynomial of degree 4 in U and of degree 2 in V. Hence we 
expect to be able to solve dyP(U, V) = 0. (In fact dyP is equal to V 
times irreducible quadratic polynomial in U.) Therefore we have exactly 
one solution V = 0, or N = M. But this is a boundary of the region, so we 
only have to check the boundary values. 

This time the boundary conditions are: M = N, M = (N + l)/2 and 
N = 2. After changing variables to U and V these become V = 0, U = 3V 
and U + V = 1, respectively. Each time we get a polynomial of degree 4. 
Thus we proceed as in the previous section. 

• P(U,0) = for U = (double root), U « 5.65 and U w -0.24; 
P(0,0) = and P(1,0) < 0, 

• P(3V,V) = for V = (double root), V « 0.55 and V ps -0.04; 
P(0, 0) = and P(3/4, 1/4) < 0, 

• P(l - V,V) = for V ss -0.52 and V w 0.29 (2 complex roots); 
P(l, 0) < and P(3/4, 1/4) < 0. 

Hence the inequality is true on the boundary and so also inside of the region. 

5. Case: 1 < N < 2 and 1 < M < (N + l)/2 

Here we take e = 1, (3 = and a = j = (A r +M-2)/ v / 2. Even though the 
symmetry described in the previous section does not exist here, we will still 
use the same rotated coordinates U = (M + N)/2 - 1 and V = (N - M)/2. 
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This time the inequality (|2.9j) becomes: 
(5.1) 

>P(U, V) = 32133332[/ 2 (-1898955433- 549628092V6 + 103783680V2tt 
- 22702680 V3tt + 345945600tt 2 - 1063944882J7 + 222614730V6E7 
+ 259459200V2tt?7 - 136216080V37r?7 + 115315200tt 2 ;7 - 531972441C/ 2 
+ 113513400v / 37rt/ 2 + 172972800tt 2 L/ 2 ) 

- 64266664(824442138V6 - 155675520%/2tt - 2201993543/7 

+ 158918760V37r[/ + 403603200tt 2 [/) (U + U 2 ) V 
+ 3759599844 (i478400tt 2 + 10405746\/6t/ + 5765760V2tt/7 - 4546773L/ 2 

+ 4074840V3ttC/ 2 + 3449600tt 2 [/ 2 ) U 2 

Note that this is still a polynomial of degree 2 in V and therefore we 
proceed as in the previous section. Unfortunately this time the only solution 
of dyP = is V equal to a rational function of U with an irreducible 
denominator of degree 2. Hence, by plugging this into djjP = we get a 
rational equation with squared irreducible polynomial of degree 2 in the 
denominator. Thus, this equation is equivalent to the numerator being 
0. Fortunately the numerator is a solvable polynomial of degree 7 with 
4 imaginary roots and 3 real roots (0, ~ —0.18 and ~ 1.8). 

Here we have the following bounds U = 3V (equivalent to M = (iV+1) /2), 
U + V = 1 (N = 2), and U = V (M = 1). So this is a triangle with vertices 
(0, 0), (3/4, 1/4) and (1/2, 1/2). Hence neither critical point is inside of this 
region. Therefore we have to check the boundary values and like before 
this means we have to find the roots of certain polynomials of degree 4 and 
values at the endpoints. 

• P(V,V) = for V = (double root), V « -0.27 and V » 0.64; 
P(0, 0) = and P(l/2, 1/2) < 0, 

• P(3V,V) = for V = (double root), V « -0.06 and V « 0.51; 
P(0, 0) = and P(3/4, 1/4) < 0, 

• P(U,1 - U) = for U w 0.48 and U ps 0.79 (2 complex roots); 
P(l/2, 1/2) < and P(3/4, 1/4) < 0. 

Thus inequality is true. 

6. Case: M > 15 

For this case we will use a method different than all other cases. Since 
we are dealing with the triangles for which two sides are long and almost 
equal, we will estimate the eigenvalue by the eigenvalue of a circular sector 
contained in the triangle T. 

Let us denote the angle between the sides of length TV and M by 7. First 
we take the isosceles triangle with angle 7 between the sides of length M. 
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We can certainly put this triangle inside the triangle T. Since the shortest 
side of this isosceles triangle has length no larger than 1, the altitude h 
satisfies h > \JM 2 - 1/4. 

Let us denote a circular sector with angle a and radius r by S(a,r). It 
is known (see (7j), that the first eigenvalue of the sector S(a,r) is j^/ a r ~ 2 , 
where j u is the first zero of the Bessel function J u (x) of order v. 

It is clear that we can put a sector 5(7, h) inside the triangle T. Hence, 
by domain monotonicity we have 

f 

(6-1) A T < As (7ih ) = -j^. 

We need to prove that 



2 2 T 2 



(6.2) Vy < vr L 



h 2 ~ 9A 2 T ' 

We have L T = 1 + M + N > 2N and A T = sin(j)NM/2 < jNM/2. 
Therefore, it is enough to prove that 

9jl h (lNM/2) 2 

(6 " 3) (M 2 - 1/4) 2 (2iV) 2 " ' 

or that 

, , 9j 2 7 2 M 2 

(6.4) < 1. 

1 ' 16tt 2 (M 2 - 1/4) 2 " 

To find the bound for j u we will use the estimate obtained in jH] 

(6.5) j v <u- ^1/3 + ^^^-1/3 
y ' 4/2 20 

where ai ~ —2.338 is the first negative zero of the Airy function. Hence we 
have 

(6.6) ^ < 1 + 2v- 2 ^ + 2v~^. 

v 

Therefore 

gj^fM 2 / , 7 x2/3 / 7 \4/3\ 2 9M 2 



16(M 2 - 1/4)' 



( , 7) _-^,( 1 + 2(3 - +2( i)-)' 

This last expression is increasing with 7 as can be easily verified by dif- 
ferentiating. Given M, the angle 7 is maximized for the isosceles triangle, 
hence 7 < 2 sin _1 (l/2M). In order to arrive at (pOjl . it is enough to show 
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It is easy to check that the function on the left side is decreasing with M, 
and that for M = 15 inequality is true. Hence this is true for any triangle 
with M > 15. 

Note also that if M — > oo, then the whole expression tends to 9/16. This 
shows that the constant 16 in the lower bound in Theorem II .11 is optimal. 

7. Script in Mathematica 

Here we give the script written in Mathematica to handle all the cumber- 
some calculations included in Sections 2 to 5. It is important to note that 
all the calculations are done symbolically. Only the exact values of the roots 
of all the polynomials are at end converted to numerical form. 

(* Section 2 *) 

(* isosceles triangle with vertices (0,0), (1,0) and (Sqrt[3],0) *) 
g[x_,y_]=Sin [Sqrt [3] \ [Pi] y]Sin[\[Pi] x/3] + \ 

Sin [\ [Pi] y/Sqrt[3]]Sin[5\[Pi] x/3] + \ 

Sin[2\[Pi] y/Sqrt[3]]Sin[4\[Pi] x/3]; 
(* other right triangles *) 
g2[x_,y_]=g[l-x,y] ; 
g3 [x_ , y_] =g [Sqrt [3] y , Sqrt [3] x] ; 

(* test functions obtained from right triangles *) 

\ [CurlyPhi] l=g [x- (a y /b) , Sqrt [3] y/b] ; 

\[CurlyPhi]2=g2[x-((a-l) y /b) , Sqrt [3] y/b] ; 

\ [CurlyPhi] 3=g3 [x- (a y /b) , y/ (Sqrt [3] b) ] ; 

(* equilateral triangle after linear transformation *) 

\[Phi] :=Sin[2\[Pi]y/b]-Sin[2\[Pi] (x+(l-a)y/b)] +Sin [2\ [Pi] (x-a y/b)] ; 

(* final test function *) 

\[Psi]=\ [Alpha] \ [CurlyPhi] 1 + \[Beta] \ [CurlyPhi] 2 + \ 

\ [Gamma] \ [CurlyPhi] 3 + \[Epsilon] \ [Phi] ; 
grad=Simplify [Integrate [D[\ [Psi] ,x] ~2+D [\ [Psi] ,y] ~2,{y,0,b>, \ 

{x,a y/b, (a-1) y/b+1}]] ; 
int=Simplify [ Integrate [\ [Psi] ~2,{y,0,b}, {x, a y/b , (a-l)y/b +1}]]; 
(* we have to prove that this is <= *) 

in=9b~2grad-4\ [Pi] ~2(1+Sqrt [a~2+b~2] +Sqrt [(a-1) ~2+b~2] ) ~2int ; 
(* change from (a, b) to (M, N) and cancel b *) 

in2=Simplify[in/b /. b~2 -> M~2 - a~2 /. a -> (M~2 - N~2 + l)/2, \ 

(N > 0) && (M > 0)]; 
(* inequality (2.9) *) 
Simplify [308788467187200in/b] 

(* Section 3 *) 

W=in2/. \[Epsilon] -> /. \ [Gamma] -> -1/6 /. \ [Beta] -> /. \ 

\ [Alpha] -> 1; 
(* Inequality (3.1) *) 
Apart [1383782400W] 
(* Critical point *) 

Reduce [(D[W, M] == 0) && (D[W, N] == 0), {M, N}] // N 
(* Boundary : roots and endpoints *) 
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Reduce [W == / . N -> 2] // N 
Reduce [W == / . M -> N - 1] // N 
Reduce [W == / . M -> N] // N 
Reduce [W == / . M -> 15] // N 
W /. M -> {1, 2} /. N -> 2 // N 
W / . M -> 15 /. N -> {15, 16} // N 



(* Section 4 *) 

W=in2/. \[Epsilon] -> 1 /. \ [Gamma] -> /. \[Beta] -> \ [Alpha] /.\ 

\ [Alpha] -> (N + M - 2)/2; 
pol = W/. M -> U - V /. N->U+V/. U->U+1; 
(* inequality (4.1) *) 
Apart [22056319084800pol , V] 
(* Critical point *) 
Reduce [D [pol , V] ==0, V] // N 
(* Boundary : roots and endpoints *) 
Reduce [pol == / . V -> 0] // N 
Reduce [pol == /. U -> 1 - V] // N 
Reduce [pol == / . U -> 3V] // N 
pol /. V -> /. U -> {0, 1} // N 
pol /. V -> 1/4 /. U -> 3/4 // N 

(* Section 5 *) 

W=in2/. \[Epsilon] -> 1 /. \ [Beta] -> /. \ [Gamma] -> \ [Alpha] /.\ 

\ [Alpha] -> (M + N - 2)/Sqrt[2]; 
pol = W /. M -> U - V /. N->U+V/. U->U+ 1; 
(* inequality (5.1) *) 
Apart [9609600pol, V] 
(* Critical points *) 
Vs = Solve [D [pol, V] ==0, V]; 
Reduce [D [pol, V] ==0, V, Reals] 
(* denominator with complx roots only*) 

Reduce [Denominator [Together [D [pol, U] /. Vs]] == 0] // N 
(* polynomial of degree 7 in U *) 
Reduce [Numerator [Together [ 

D[pol, U] /. Vs]] == 0] // N 
(* Boundary : roots and endpoints *) 
Reduce [pol == / . U -> 3V] // N 
Reduce [pol == / . U -> V] // N 
Reduce [pol == / . V -> 1 - U] // N 
pol /. U -> 1 - V /. V -> {1/4, 1/2} // N 
pol /. U -> /. V -> // N 
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